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ABSTRACT 

We report on the long term X-ray monitoring of the outburst decay of the low magnetic field mag- 
netar SGR 0418+5729 , using all the available X-ray data obtained with RXTE, Swift, Chandra, and 
XMM-Newton observations, from the discovery of the source in June 2009, up to August 2012. The 
timing analysis allowed us to obtain the first measurement of the period derivative of SGR 0418+5729: 
P = 4(1) x 10~ 15 ss" 1 , significant at ~ 3.5a confidence level. This leads to a surface dipolar magnetic 
field of £?dip ~ 6 x 10 12 Gauss. This measurement confirms SGR 0418+5729 as the lowest magnetic 
field magnetar. Following the flux and spectral evolution from the beginning of the outburst up to 
~1200days, we observe a gradual cooling of the tiny hot spot responsible for the X-ray emission, from 
a temperature of ~ 0.9 to 0.3 keV. Simultaneously, the X-ray flux decreased by about 3 orders of 
magnitude: from about 1.4 x 10~ n to 1.2 x 10" 14 ergs -1 cm~ 2 . Deep radio, millimeter, optical and 
gamma-ray observations did not detect the source counterpart, implying stringent limits on its multi- 
band emission, as well as constraints on the presence of a fossil disk. By modeling the magneto-thermal 
secular evolution of SGR 0418+5729 , we infer a realistic age of ^550 kyr, and a dipolar magnetic field 
at birth of <~ 10 14 G. The outburst characteristics suggest the presence of a thin twisted bundle with 
a small heated spot at its base. The bundle untwisted in the first few months following the outburst, 
while the hot spot decreases in temperature and size. We estimate the outburst rate of low magnetic 
field magnetars to be about one per year per galaxy, and we briefly discuss the consequences of such 
result in several other astrophysical contexts. 

Subject headings: sources (individual): SGR 0418+5729 — stars: magnetic fields — stars: neutron 



1. INTRODUCTION 

Neutron stars showing magnetar-like activity (compris- 
ing the anomalous X-ray pulsars, soft gamma repeaters 
and a high magnetic field pulsar) are a small group of 

1 Institute of Space Sciences (CSIC-IEEC), Campus UAB, 
Faculty of Science, Torre C5-parell, E-08193 Barcelona, Spain 

2 INAF/Osscrvatorio Astronomico di Roma, via Frascati 33, 
1-00040 Monteporzio Catone, Italy 

3 Department de Fisica Aplicada, Universitat dAlacant, Ap. 
Correus 99, E-03080 Alacant, Spain 

4 Universita di Padova, Dipartimcnto di Fisica e Astronomia, 
via F. Marzolo 8, 1-35131 Padova, Italy 

5 Mullard Space Science Laboratory, University College Lon- 
don, Holmbury St. Mary, Dorking, Surrey RH5 6NT, UK 

6 INAF/Istituto di Astrofisica Spazialc c Fisica Cosmica - Mi- 
lano, via E. Bassini 15, 1-20133 Milano, Italy 

7 JILA and Department of Astrophysical and Planetary Sci- 
ences, University of Colorado, Boulder, CO 80309, USA 

8 IUSS - Istituto Universitario di Studi Superiori, Piazza della 
Vittoria 15, 1-27100 Pavia, Italy 

9 INFN Istituto Nazionale di Fisica Nucleare, Sezione di 
Pavia, Via Bassi 6, 1-27100 Pavia, Italy 

10 Universita dcgli Studi di Pavia, Dipartimcnto di Fisica Nu- 
cleare e Teorica, via Bassi 6, 1-27100 Pavia, Italy 

11 INAF/Osscrvatorio Astronomico di Cagliari, localita Pog- 
gio dci Pini, strada 54, 1-09012 Capoterra, Italy 

12 Sabanci University, Orhanh-Tuzla, 34956 Istanbul, Turkey 

13 NASA Marshall Space Flight Center, Huntsville, AL 35812, 
USA 

14 AIM (UMR 7158 CEA/DSM-CNRS-Universite Paris 
Diderot) Irfu/Scrvicc d'Astrophysique, Saclay, FR-91191 Gif- 
sur-Yvette Cedex, France 

15 Kepler Institute of Astronomy, University of Zielona Gora, 
Lubuska 2, 65-265, Zielona Gora, Poland 

16 SRON-Nethcrlands Institute for Space Research, Sorbon- 
nelaan 2, 3584 CA Utrecht, The Netherlands 



X-ray pulsars (about twenty objects) with spin periods 
between 0.3-12 s, whose strong persistent and/or flar- 
ing emission are hard to explain by the common sce- 
narios for rotation powered pulsars or accreting pulsars. 
In fact, the very strong X-ray emission of these objects 
(Lx ~ 10 35 ergs _1 ) is too high and/or variable to be fed 
by the rotational energy alone (as in the radio pulsars), 
and no evidence for a companion star has been found, 
hence ruling out accretion in a binary. Accretion from 
a fossil disk remnant of the supernova explosion might 
be responsible for part of the observational properties of 
these objects, but it fails to explain some of their char- 
acteristics, such as the flaring X-ray activity. Their in- 
ferred magnetic fields, under the typical assumption of 
magnetic dipolar losses alone, appear to be as high as 

B dip ~ 3.2 x 10 19 \/P? ~ 10 14 - 10 15 G (see Mereghetti 
2008 for a review). These strong fields are believed to 
form cither via a dynamo action in a rapidly rotating 
proto-neutron star (<3ms; Thompson & Duncan 1995), 
or they are thought to be fossil fields remnant of a highly 
magnetic massive star (<~1 kG; Ferrario & Wickramas- 
inghe 2006). Because of these high B fields, the emis- 
sion of "magnetars" is thought to be powered by the de- 
cay and the instability of their strong fields (Duncan & 
Thompson 1992; Thompson & Duncan 1993, Thompson, 
Lyutikov & Kulkarni 2002). Their powerful X-ray out- 
put is usually well modeled by thermal emission from the 
neutron star hot surface, reprocessed in a twisted magnc- 
tosphcrc through resonant cyclotron scattering (Thomp- 
son, Lyutikov & Kulkarni 2002; Nobili, Turolla & Zane 
2008; Rea et al. 2008; Zane et al. 2009), a process 
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TABLE 1 

Journal of all the X-ray observations of SGR 0418+5729. 



Instrument 




Obs. ID 


Starting 


; date 


Exp. (ks) 


Counts s 1 


Flux c 


kT BB (keV) d 


BB norm. 6 


RXTE/PCA" 


94048 


2009 06-11/11-24 


194.2 




— 


— 


- 




Swift/XRT 




00031422001 


2009-07-08 


20:48:01 


2.9 


0.229±0.008 


13.4+1.0 


0.88+0.05 


2.2±0. 


.2 


Swift/XRT (PC) 


00031422002 


2009-07-09 


00:04:01 


10.6 


0.245±0.003 


13.8±0.7 


0.94±0.03 


1.7+0 


.1 


Swift/XRT (PC) 


00031422003 


2009-07-10 


00:15:01 


5.6 


0.179+0.005 


11.0+0.7 


0.95±0.05 


1.33+0 


.13 


Swift/XRT (WT) 


00031422004 


2009-07-12 


00:27:01 


7.1 


0.218+0.006 


11.5+0.7 


0.91+0.04 


1.68±0. 


.14 


Chandra/HRC-P 


10168 


2009-07-12 


06:06:43 


24.1 


0.317+0.005 






- 




Swift/XRT (WT) 


00031422006 


2009-07-15 


00:48:39 


7.7 


0.252±0.006 


■ 


0.93±0.03 


1.69±0. 


.12 


Swift/XRT (WT) 


00031422007 


2009-07-16 


00:53:01 


16.4 


0.217±0.004 


11.7+0.8 


0.93+0.02 


1.57+0 


.08 


XMM-Newton/EPIC* 


0610000601 


2009-08-12 


21:09:12 


67.1 


1.281±0.005 


6.75+0.07 


0.897±0.007 


1.038+0 


.018 


Swift/XRT 


(PC) 


00031422008 


2009-09-20 


21:09:00 


9.4 


0.066+0.002 


3.66±0.30 


0.82±0.05 


0.79±0. 


.09 


Swift/XRT 


(PC) 


00031422009 


2009-09-22 


00:43:00 


7.6 


0.072±0.003 


3.58±0.40 


0.82±0.05 


0.79+0 


.11 


Swift/XRT 


(PC) 


00031422010 


2009-11-08 


00:36:01 


15.1 


0.043±0.002 


2.14+ 0.20 


0.82±0.05 


0.47+0 


.05 


Swift/XRT 


(PC) 


00031422011" 


2010-01-14 


08:06:01 


3.6 


o.oig+o.ooi 6 


1.05+0.10 


0.75+0.07 


0.32+0 


.06 


Swift/XRT 


(PC) 


00031422012 61 


2010-01-15 


13:08:01 


3.7 


o.oig+o.ooi 6 






" 




Swift/XRT 


(PC) 


00031422013 61 


2010-01-16 


08:14:01 


4.0 


0.019±0.001 b 






" 




Swift/XRT 


(PC) 


00031422014 61 


2010-01-17 


06:47:01 


3.8 


0.019±0.001 b 










Swift/XRT 


(PC) 


00031422015 62 


2010-02-14 


17:33:01 


4.5 


0.0172±0.0008 c 


0.76+0.11 


0.74+0.04 


0.25+0 


.03 


Swift/XRT 


(PC) 


00031422016 62 


2010-02-15 


17:37:01 


4.5 


0.0172±0.0008 c 










Swift/XRT 


(PC) 


00031422017 M 


2010-02-16 


01:38:01 


4.6 


0.0172±0.0008 c 










Swift/XRT 


(PC) 


00031422018 62 


2010-02-17 


09:49:01 


1.6 


0.0172±0.0008 c 




" 


" 




Swift/XRT 


(PC) 


00031422019 62 


2010-02-18 


16:14:01 


3.9 


0.0172±0.0008 c 










Swift/XRT 


(PC) 


00031422020 62 


2010-02-19 


00:23:01 


3.2 


0.0172±0.0008 c 






" 




Swift/XRT 


(PC) 


00031422021 63 


2010-07-09 


06:50:01 


3.6 


0.0023±0.0003 d 


0.10+0.04 


0.47+0.13 


0.23±0 


.13 


Swift/XRT 


(PC) 


00031422022 63 


2010-07-10 


18:11:00 


5.2 


0.0023±0.0003 d 










Swift/XRT 


(PC) 


00031422023 63 


2010-07-11 


05:19:01 


5.0 


0.0023±0.0003 d 










Swift/XRT 


(PC) 


00031422024 63 


2010-07-11 


23:06:01 


5.4 


0.0023±0.0003 d 










Swift/XRT 


(PC) 


00031422025 63 


2010-07-13 


00:47:01 


4.9 


0.0023+0. 0003 d 










Chandra/ACIS-S 


12312 


2010-07-23 


15:04:09 


30.0 


0.0017±0.0008 


0.13+0.02 


0.68±0.04 


0.061+0 


.008 


XMM-Newton/EPIC 


0605852201 


2010-09-24 


01:54:56 


34.2 


0.0370±0.0020 


0.16±0.02 


0.69±0.05 


0.07±0. 


.01 


Chandra/ACIS-S 


13148 


2010-11-29 


05:59:57 


30.0 


0.0038+0.0004 


0.021+0.002 


0.38±0.11 


0.12+0 


.06 


XMM-Newton/EPIC 


0672670201 


2011-03-10 


03:15:53 


35.0 


0.0071±0.0007 


0.015+0.002 


0.32±0.05 


0.21+0 


.08 


Chandra/ACIS-S 


13235 


2011-07-20 


02:26:12 


77.0 


0.0033±0.0002 


0.015+0.003 


0.37±0.04 


0.11+0 


.02 


XMM-Newton/EPIC M 


0672670401 


2011-09-09 


15:27:23 


33.0 


0.0071+0.0006 6 


0.016+0.002 


0.28±0.05 


0.34+0 


.13 


XMM-Newton/EPIC M 


0672670501 


2011-09-11 


21:47:41 


48.5 


0.0071+0.0006 6 










Chandra/ACIS-S 


13236 


2011-11-26 


11:48:02 


75.0 


0.0026+0.0002 


0.015+0.002 


0.35+0.07 


0.13+0 


.05 


XMM-Newton/EPIC* 


0693100101 


2012-08-25 


14:18:08 


78.2 


0.0058+0.0004 


0.012+0.001 


0.32+0.05 


0.16+0 


.05 



a The RXTE-PCA and the Chandra HRC-I were used only for the timing analysis. 

b These observations were merged in the timing and spectral analysis to improve statistics. 

c Absorbed flux in the 0.5-10 keV energy range, and in units of 10~ 12 erg s _1 cm — 2 . Errors in the table are at 90% confidence level. 

d Fitted model is: phabs*bbodyrad; N H = (1.15 ± 0.06) X 10 21 cm" 2 and \l = 1-19 (for 940 dof). 

e The BB radius in km is the square root of this BB normalization, times the distance in units of 10 kpc . 

* See 33.1|for details on the modeling of these observations. 



favored under these extreme magnetic conditions. On 
top of their persistent X-ray emission, magnetars emit 
very peculiar flares and outbursts on several timescales, 
from fractions of a second to years reaching very high, 
super-Eddington luminosities (10 38 — 10 46 ergs _1 ). These 
flares are most probably caused by rearrangements of the 
twisted magnetic field lines, either accompanied or trig- 
gered by fractures of the neutron-star crust (Thompson 
& Duncan 1995; Perna & Pons 2011). 

Transient events are a characteristic signature of mag- 
netar emission, and one of the main ways to discover new 
sources of this class and study their physics. From the 
discovery of the first transient less than a decade ago, we 
now count about a dozen of outbursts, which increased 
the number of known magnetars by a third in six years 
(see Rea & Esposito 2011; Rea 2013 for recent reviews). 
Magnetar outbursts might involve their multi-band emis- 
sion resulting in an increased activity from radio to hard 



X-ray, usually with a soft X-ray flux increase of a factor 
of 10-1000 with respect to the quiescent level. An as- 
sociated X-ray spectral evolution is often observed, with 
a spectral softening during the outburst decay (Rea et 
al. 2009). The flux decay timescale varies substantially 
from source to source, ranging from a few weeks to sev- 
eral years (Rea & Esposito 2011; Pons & Rea 2012). 

The extensive follow-up of magnetars undergoing an 
outburst yielded the most unexpected discovery of the 
past years in the magnetar field. Prompted by the de- 
tection of typical magnetar-like bursts and a powerful 
outburst, a new transient magnetar with a spin period 
of -9 s was discovered in 2009, namely SGR 0418+5729 
(van der Horst et al. 2010, Esposito et al. 2010). How- 
ever, after more than 2 years of extensive monitoring, no 
period derivative was detected. This led to an upper limit 
on the source surface dipolar field of Bdi P < 7.5 x 10 12 G 
(Rea et al. 2010). For the first time, we detected a mag- 
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Fig. 1. — Left panel: Evolution of the pulse phases with time (upper panel). The solid lines represent the timing solution without (linear) 
and with a P component (quadratic). The time residuals (lower panel) are relative to the quadratic fit. Right panel: Pulsed fraction 
evolution in the 0.5— lOkeV band for the Swift (red), Chandra (green) and XMM-Newton (blue) observations. 



netar with a low dipolar magnetic field, showing that a 
critical magnetic field is not necessary for a neutron star 
in order to display magnetar-like activity. In turn, this 
means that many seemingly normal pulsars could turn 
out as magnetars at anytime (this was supported by the 
discovery of a second low-B magnetar followed soon af- 
ter; Rea et al. 2012; Scholz et al. 2012). After the 
discovery of this low dipolar magnetic field soft gamma 
repeater, several models were put forward to explain its 
puzzling emission. They involve the possible presence of 
a fall-back disk slowing down the pulsar up to the cur- 
rent spin period (Alpar et al. 2011), a tiny inclination 
angle between the magnetic and rotational axis resulting 
in a higher inferred magnetic field (Tong & Xu 2012), a 
pulsar with a strongly magnetized core (Soni 2012), an 
old quark nova (Ouyed, Leahy & Niebergal 2011), or a 
massive highly magnetized, rotating white dwarf (Mal- 
heiro, Rueda & Ruffini 2012). In Rea et al. (2010) and 
Turolla et al. (2011), we suggested that a non-dipolar 
component of the field, larger than the measured dipolar 
one, can be responsible for the behavior of this magnetar, 
if it has a relatively old age (~1 Myr). 

In this paper we present the complete study of the 
outburst of the low dipolar magnetic field magnetar 
SGR 0418+5729, from the first outburst phases until 
about 3 years after its onset. This long term monitor- 
ing campaign using several X-ray satellites, allowed us to 
estimate SGR0418+5729's period derivative, and follow 
the cooling of its surface temperature during the outburst 
decay up to the (probable) quiescent level. Furthermore, 
we inferred limits on its emission in the radio, millimeter, 
optical and gamma-ray bands. We discuss our findings in 
terms of the magneto-thermal history of this magnetar, 



discuss the current limits on the presence of a fossil disk, 
and present some discussion on the broader consequences 
of the discovery of low magnetic field magnetars. 

2. X-RAY OBSERVATIONS AND DATA REDUCTION 

In this study, we used data obtained from several differ- 
ent satellites (see Table 1 for a summary). We describe 
below the observations and data analysis. Part of the 
data we used in this paper were already published by 
van der Horst et al. (2010), Esposito et al. (2010) and 
Rea et al. (2010). 



2.1. 



Swift data 

The X-Ray Telescope (XRT; |Burrows et aT|[2005| ) on- 



board Swift uses a front-illuminated CCD detector sensi- 
tive to photons between 0.2 and 10 kcV. Two main read- 
out modes are available: photon counting (PC) and win- 
dowed timing (WT). PC mode provides two dimensional 
imaging information and a 2.5073 s time resolution; in 
WT mode only one-dimensional imaging is preserved, 
achieving a time resolution of 1.766 ms. The XRT data 
were uniformly processed with xrtpipeline (version 12, 
in the heasoft software package version 6.11), filtered 
and screened with standard criteria, correcting for effec- 
tive area, dead columns, etc. The source counts were 
extracted within a 20-pixel radius (one XRT pixel cor- 
responds to about 2'.'36). For the spectroscopy, we used 
the spectral redistribution matrices in CALDB (20091130; 
matrices version v013 and v014 for the PC and WT data, 
respectively), while the ancillary response files were gen- 
erated with XRTMKARF, and they account for different 
extraction regions, vignetting and point-spread function 
corrections. 
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Fig. 2. — Spectral modeling of all observations listed in Table 1. Left and middle panels: spectra and residuals for all observations ( Swift 
(red), Chandra (green) and XMM-Newton (blue))fitted simultaneously with a single blackbody model (left) and using an RCS model only 
for the first XMM-Newton observation (center). Right panel: unfolded spectrum relative to the modeling shown in the central panel. 



2.2. RXTE data 



The Pro portional Counter Array (PCA; jJahoda 
1996) on-board RXTE consists of five collimated 



et al 



xenon/ methane multi-anode Proportional Counter Units 
(PCUs) operating in the 2-60 keV energy range. Raw 
data were reduced using the FTOOLS package (version 
6.11). To study the timing properties of SCR 0418+5729, 
we restricted our analysis to the data in Good Xenon 
mode, with a time resolution of 1 /Lis and 256 energy bins. 
The event-mode data were extracted in the 2-10 keV en- 
ergy range from all active detectors (in a given observa- 
tion) and all layers, and binned into light curves of 0.1s 
resolution. We use here 46 RXTE/PCA observations of 
SGR 0418+5729, spannig the first 6 months of the out- 
burst, until the source flux decayed below the instrument 
detection level. The total 194.2 ks exposure time is di- 
vided in observations of 0.6 to 13.6 ks exposure each. See 
Esposito et al. (2010) for further details on the Swift and 
RXTE observations. 

2.3. Chandra data 

The Chandra X-ray Observatory monitored 
SGR 0418+5729 five times during the past three 
years. The first one with the High Resolution Imaging 
Camera (HRC-I; Zombeck et al. 1995) and the following 
four observations with the Advanced CCD Imaging 
Spectrometer (ACIS-S; Garmire et al. 2003). Data 
were analyzed using standard cleaning procedure^ | and 
CIA0 version 4.4. The HRC-I camera does notnave 
a sufficient spectral resolution, and it was used only 
for the timing analysis; it has a timing resolution of 
^16 /xs. All ACIS-S observations were performed in 
VERY FAINT mode, with only the S7 CCD on, resulting 
in a timing resolution of 0.44 s. Photons were extracted 
from a circular region with a radius of 3" around the 
source position, including more than 90% of the source 
photons, and background was extracted from a similar 
region far from the source position. 

2.4. XMM-Newton data 
SGR 0418+5729 was observed six times with XMM- 



Newton ( Jansen et al.||200T ). Data have been processed 

17 http://asc.harvard.edu/ciao/threads/index.html 



using SAS version 12, and we have employed the most 
updated calibration files available at the time the re- 
duction was performed (August 2012). Standard data 
screening criteria are applied in the extraction of scien- 
tific products. For our spectral analysis we used only the 
EPIC-pn camera (Turner et al. 2001) which provides the 
spectra with the best statistics, while the MOS cameras 
(Striider et al. 2001) were added in the timing analysis. 
The EPIC-pn camera was set in Small Window (timing 
resolution of 6 ms) and Full Frame (73 ms) modes in the 
first two observations, respectively, and in Large Window 
mode for all the following ones (48 ms) , with the source 
at the aim-point of the camera, and the MOS cameras 
in Small Window mode (0.3 s). We extracted the source 
photons from a circular region of 30"radius, and a similar 
region was chosen for the background in the same CCD. 
We restricted our spectral analysis to photons having 
PATTERN < 4 and FLAG=0 for the EPIC-pn data. 

3. RESULTS OF THE X-RAY MONITORING 
3.1. X-ray spectral modeling 

For the spectral analysis we used source and back- 
ground photons from the Swift, Chandra and XMM- 
Newton observations extracted as described in the pre- 
vious section (we also checked our results using a larger 
extraction region for the background spectra). The re- 
sponse matrices were built using ad-hoc bad-pixel files 
built for each observation. We used the XSPEC package 
(version 12.4) for all fittings, and the phabs absorption 
model with the Anders & Grevesse (1989) abundances, 
and the Balucinska-Church & McCammon (1992) photo- 
electric cross-sections. We restricted our spectral model- 
ing to the 0.7-10 keV energy band, excluding bad chan- 
nels when needed. The Swift spectra were binned in 
order to have at least 30 counts per spectral bin. XMM- 
Newton spectra were grouped such to have at least 100, 
50 and 40 photons per bin in the first three observations, 
respectively, and a minimum of 30 counts in the sub- 
sequent observations. On the other hand, all Chandra 
spectra have at least 25 counts per bin. 

We started the spectral analysis by fitting all the 
spectra together (see Table [I]) with a single component 
model: an absorbed blackbody or a powerlaw model. 
While the former gave a good fit, a single powerlaw could 
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Fig. 3. — Spectral evolution with time. Top panel: flux evolution 
for the absorbed 0.5-10 keV flux (black), and for the bolometric 
unabsorbed flux (red). Middle and bottom panels: evolution of 
the blackbody temperature and radius, calculated at infinity (the 
latter assuming a 2kpc distance). 




Fig. 4. — Fitted blackbody models (see Table 1), with the first 
and last observations labeled as red squares (see Figure 3 and text 
for details). 



not reproduce all the spectra at the same time. Fixing 
the absorption value to be the same for all spectra, for 
a single blackbody model (phabs*bbodyrad) we find an 
acceptable fit with iV H = (1.15 ± 0.06) x 10 21 cm" 2 and 
xl = 1.19 (940 dof; errors on the spectral parameters are 
all reported at 90% confidence level). However, not un- 
expectedly, the best collected spectrum (the first XMM- 
Newton observation on 2009-08-12; see Table 1) gave 
bad residuals at lower and higher energies (see Figure 
[2] left panel). We tried to model this observation alone, 
and indeed a single absorbed blackbody or power-law 
components were not reproducing this spectrum prop- 
erly (xl > 2). We then used a composite model. Good 
fits were found both using an absorbed blackbody plus 
a powerlaw (phabs* (bbodyrad + power); Nn= (6.32 ± 
0.04) x 10 21 cm- 2 , kT= 0.91 ± 0.07 keV, T = 2.82 ± 0.16, 
and xl = 0.97 for 392 dof), and an absorbed reso- 
nant cyclotron scattering model (RCS: Rea et al. (2007, 
2008), or NTZ: Zane et al. (2009)). Two blackbodies 
were also producing acceptable reduced chi-square val- 
ues (xl = 1.01 for 392 dof) but with worse residuals 
at higher energies (this is compatible with what found 
in Esposito et al. (2011) and Turolla et al. (2011)). 
The parameters we found for the resonant cyclotron 
scattering models are: Nu~ (1.9 ± 0.3) x 10 21 cm~ 2 , 
r = 8.8 ± 1.2, p = 0.21 ± 0.08 and kT= 0.63 ± 0.11 keV 
(xl = 1-08 for 392 dof) for the RCS model; and N K = 
(1.8±0.2) x 10 21 cm- 2 , Acf> = 1.9±1.0, ^ bu ik = 0.13±0.05 

and kT= 0.88 ± 0.07 keV (xl = L11 for 392 dof ) for the 
NTZ model. 

We then continued our spectral modeling of all spec- 
tra together by adding a further component only for this 
observation (adding a further component for all spec- 
tra was not significantly changing the goodness of the 
fit; xl = 1-13 (938 dof); see Figure^ middle and right 
panels). In Table [T] we report the values of the single 
absorbed blackbody model (see also Figure[2] left panel, 
and Figure|4| , since when fitting all the spectra by using a 
composite model only for the first XMM-Newton obser- 
vation, we find no change in the parameters of the other 
spectra with respect to the single blackbody fit. How- 
ever, although a blackbody plus powerlaw model gives 
a good fit when fitting the first XMM-Newton observa- 
tion alone, it is not so when fitting all data together. 
This is because the powerlaw component produces an 
unrealistic N-r increase, which does not match the value 
required by all the other observations modeled by a sin- 
gle blackbody. We then use one of the resonant cyclotron 
scattering models, the RCS model, for the joint-fit as an 
empiric model for the first XMM-Newton observation pj 

In addition to the joint-spectral modeling, we also nt- 
ted all the spectra individually. Beside the first XMM- 
Newton observation discussed above, the last XMM- 
Newton observation, when fitted alone with a single 
blackbody model, did not give a good chi-square (xl = 
2.2 for 16 dof). Given the low number of counts in 



18 Note that both the RCS and NTZ models are built for higher 
surface dipolar fields, hence the fact that they provide a very good 
fit to the data is probably just an indication of the presence of some 
magnetospheric distorsion. However, no real physical information 
can be derived from the resulting magnetospheric parameters. For 
the purpose of this work, we are mainly interested in the surface 
thermal cooling of the source. 
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Fig. 5. — Pulse profiles evolution in the 2-10 keV (RXTE) and 0.5— lOkeV (Chandra and XMM-Newton) energy ranges, for most of the 
observations reported in Table 1. Epoch increases from left to right, and top to bottom. 

vation (see Table 1) and leave the second blackbody free 
to vary. In neither of those two cases we could find any 
improvement in the modeling of the data. Note that al- 
though a joint two blackbody model can fit the first few 
observations (Turolla et al. 2011), this is no longer the 
case when modeling together all the data collected in the 
whole 1200 days long timespan. 

3.2. X-ray timing analysis 

All the Chandra and XMM-Newton event files col- 
lected between November 2010 and August 2012 were 
used in order to extend the coherent timing solution 
we derived in Rea et al. (2010): P = 9.07838827(4) s 
90% el., and 3c first period derivative upper limit of 



the spectrum of this observation (^400 background- 
subtracted counts), this deviation from the blackbody 
model had only a marginal effect on the joint fit. A bet- 
ter fit was found adding a second blackbody (xi = 1-2 
for 14 dof ) or (with a slightly worse chi-square) a power- 
law component (xt = 1-5 for 14 dof). However, given 
the reduced number of counts collected in this observa- 
tion, a detailed modeling of the quiescent spectrum of 
SGR 0418+5729 will be possible only when more data 
will be accumulated. 

We also tried to: 1) model all the spectra with two 
blackbodies leaving one of the blackbodies with a fixed 
area mimicking the whole surface emission, and 2) fix 
one blackbody to the value observed in the last obser- 
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TABLE 2 

Pulse phase spectroscopy of the first XMM-Newton 

OBSERVATION OF SGR 0418+5729. 



Fig. 6. — Pulse profile (normalized counts/s versus phase) as a 
function of energy, relative to the first XMM-Newton observation 
(see Table 1). 

\P\ < 6.0 x KT^ss- 1 at epoch 54993.0 MJD. Photon 
arrival times were corrected to the barycentre of the solar 
systerrp*] Timing analysis was carried out by means of 
a phase- fitting technique (deta ils on this technique are 
given in IDall Osso et al.|[2003| see also Esposito et al. 
2010 for further details on this source) . Given the intrin- 
sic variability of the pulse shape as a function of time 
(see Figurepl, we inferred the phase of the modulation 
by fitting the average pulse shape of each observation 
with a number of harmonics, the exact number of which 
is variable and determined by requesting that the inclu- 
sion of any higher harmonic is statistically significant 
(by means of an F-test) . All data reported in Table 1 
were folded using a reference period 9.07838880562798 s 
at epoch 54993 MJD, and fitted with one or more har- 
monics. In Figure]]] we plot the phases at which the fun- 
damental sine function is equal to zero in its ascending 
part (positive derivative). 

The fit of the resulting pulse phases with a linear com- 
ponent gives a reduced xl ~ 3.2 for 26 degree of freedom 
(d.o.f. hereafter). The inclusion of a quadratic term, cor- 

19 We have corrected the arrival times of the last XMM- 
Newton observation for the 2012 June 30 leap second (see 
http:/ /xmm.esa.int/sas/current/watchout/12.0.0/leapsec_2012.shtml 
tor further details J . 
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3.11 ± 0.21 



Absorbed flux in the 0.5-10 keV energy range, and in u nits 



^(lO -3 photons cm- 2 s- 1 ). See also p3 



Fitted model is: 



0.5) X 10 



21, 



phabs* (bbodyrad+power) ; 7Vh= (6.9 ± 
kT SB = 0.91 ± 0.01 keV, BB norm= 



0.82 ± 0.05, and xl = 0.96 (for 680 dof). 

responding to a first period derivative component, was 
found to be significant at a confidence level of 3.5c (by 
means of a F-test). The resulting best-fit solution corre- 
sponds to P = 9.07838822(5) s (la c.L, 2 parameters of 
interest; epoch 54993.0 MJD) and P = 4(1) x 10~ 15 ss" 1 
with a reduced Xr ~ 2.1 (for 25 d.o.f.; see also Figure]!]). 
The new timing solution implies a r.m.s. variability of 
only 0.2 s . As depicted above, the time evolution of 
the phase can be described by a relation of the form 
<p = <po + 2n(t - t )/P - tt(< - t )2P/P 2 . 

To further assess the significance of the quadratic com- 
ponent, reflecting the period derivative, we performed 
detailed Monte Carlo (MC) simulations assuming as the 
null model a simple linear relation (see Protassov et al. 
2002 for further details). By running 10 5 MC simulations 
we verified that the quadratic component is significant at 
> 99.96% confidence level, which is in very good agree- 
ment with what estimated by means of the F-test. We 
notice that examining the simulated data with different 
sampling distributions does not change significantly the 
results. Furthermore, we performed the same MC simu- 
lations using the whole set of observations but without 
considering the last XMM-Newton observation. We find 
a chance probability for the addition of a quadratic com- 
ponent of 0.65% (< 3cr). Given the strong influence of 
the last of our observations in the determination of the 
period derivative, we will perform further X-ray obser- 
vations in the next few years in order to increase the 
significance of the current P measurement. 

Using this P measurement, we infer a surface dipo- 
lar magnetic field strength of Bd ip = (6 ± 2) x 10 12 G, 
calculated at the neutron star equator. This val ue is 
fully c o nsiste nt with the 3er upper limit reported in |Rea| 



et al. 



2010[ ). We also estimate a characteristic age 
of t c ~ P/2P ~ 35Myr, and a rotational power of 
£- 3.9 x 10 46 P/P 3 - 2 x lO^ergs" 1 . 

Based on the above phase coherent timing solution, we 
also studied the pulse shape and pulsed fraction evolu- 
tion. Figure]!] shows the pulsed fraction evolution as a 
function of time. There is an evident increase starting 
soon after the burst detection with a recovery towards 
an asympthotic quiescent value which appears to be at 
about the 70-80% level. 

In Figures]!] [7] and [6] we study in the detail the shape 
of the pulse profile as it evolves in time and in energy. By 
looking at the profile shapes of all X-ray observations per- 
formed so far, the source appears to be switching among 
a three/two/single peak shape during the early outburst 
phases, with no clear trend in time (Figures]!]). The pulse 
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profile stabilizes to a single peak about three months af- 
ter the outburst onset. However, studying in detail the 
first and the last XMM-Newton observations, a few key 
pieces of information can be extracted: a) at lower ener- 
gies (<lkeV), the pulse profile is mainly single peaked, 
while the second, and possibly also the third peak ap- 
pears at higher energies (see Figure[6| ; b) the main com- 
ponent of the pulse profile continues to be at the same 
phase over the whole outburst decay (see Figure[7|. 

3.3. Pulse Phase Spectroscopy 

We performed a pulse phase spectroscopy of the first 
XMM-Newton observation. A clear pulse phase depen- 
dence of the spectrum is already observed by simply look- 
ing at the pulse profile changes as a function of the energy 
(see Figure[6]) . In order to quantify the spectral variabil- 
ity as a function of the rotational phase, we performed 
a pulse phase spectroscopy extracting the spectra from 
phases 0-0.4, 0.4-0.6, and 0.6-1. These phase intervals 
were chosen by looking at Figure|6] in order to isolate 
the dip in the 1-4 keV pulse profiles at phase ^0.55. In 
Table 2 we report the results of our modeling. The phase- 
averaged spectrum is not well fit by either a single black- 
body nor a powerlaw. We then used an absorbed black- 
body plus powerlaw modeling (note that the RCS and 
NTZ models are not suited for phase resolved analysis 
since they are intrinsically phase- average), using for the 
photoelectric absorption model the same cross-section 
and abundances as for the phase-average spectrum (see 
£3.1). The blackbody temperature and radius were con- 



sistent in all three spectra, hence we fixed them to be 
the same for all spectra (kT BB = 0.91 ± 0.01 keV, BB 
norm= 0.82 ± 0.05), while a variability > 3cr has been 
observed in the photon index (it changed from about 2.9 
to 4.1 between the spectra of the first peak and the dip). 

However, from Figure]?] it is clear that the main dif- 
ference in the spectra is at lower energies. In particu- 
lar, above 5 keV the three spectra are very similar, while 
the 0.4-0.6 phase-resolved spectrum seems to have less 
counts than the other below such energy. 

4. GREEN BANK TELESCOPE RADIO OBSERVATIONS 

We observed SGR 0418+5729 using the 101m Green 
Bank Telescope (GBT) on 2012 October 4th, during the 
return to quiescence. Data were acquired with the Green 
Bank Ultimate Pulsar Processing Instrument (GUPPI; 
DuPlain ct al. 2008) at a central frequency of 2.0 GHz 
(with a bandwidth of 800 MHz, integration time of ~ 
5400 s and sampling time of 655 fis) and 820 MHz (with 
a bandwidth of 200 MHz, integration time of ~ 5600 s, 
and sampling time of 655 (is). To minimize the dispersive 
effects of the interstellar medium, the bandwidths were 
split into 2048 and 512 channels, respectively. The work- 
ing of the system was checked looking at the pulsar PSR 
B0450+55. A mask was first applied to the full resolu- 
tion data for reducing the effects of impulsive RFI and of 
bad channels. Then the cleaned data were dowsampled a 
factor 2 in sampling time, matching the frequency resolu- 
tion in order to have a maximum dispersion smearing of 
order 1.3 ms in each channel for a pulsar with dispersion 
measure (DM) ~ 100 pc cm~ 3 

The ephemerides acquired from the X-ray observations 
(see §3.2), were used to fold the resulting data in ~ 3- 
min long subintegrations at the known magnetar period. 



We also folded the data at half, one third and a quarter 
of the nominal period in order to detect putative higher 
harmonics components of the intrinsic signal, in case the 
latter were deeply contaminated by interference (RFI). 
Folding was done using dspsr (van Straten & Bailcs 
2011). The sub-integrations and the frequency channels, 
cleaned from RFI, were then searched around the pulsar 
period P and over a wide range of DM values (from 
to 1000 pc cm~ 3 ) to find the P-DM combination maxi- 
mizing the signal-to-noise ratio. No dispersed signal was 
found in the data down to a signal-to-noise limit of 10 in 
both datasets. Given the parameters of the antenna and 
of the receiver^] and assuming a pulsar with a duty cy- 
cle of 10%, that translates in flux densities of ~ 0.02 mJy 
and ~ 0.05 mJy for the 2 GHz and 820 MHz observa- 
tions, respectively. Data were also blindly searched for 
a periodic signal in the Fourier domain, and for single 
de-dispersed pulses (within a DM range from to 200 pc 
cm -3 ). No signal was found in either the Fourier domain 
(down to a spectral signal-to-noise ratio 4) or in the sin- 
gle pulse searches (down to a signal-to-noise ratio 5 for 
the individual pulses). 

No previous search for pulsed radio emission had been 
performed at 2 GHz, whereas the flux density value at 
820 MHz improves by ~ 15% the limit of the observation 
at 820 MHz performed on 2009 July 19th (Lorimer et al. 
2009, Atel 2096), when the source was in the phase of 
outburst. Assuming a typical pulsar spectral index of 
1.7, a typical duty cycle ~ 10% and a distance of 2 kpc 
(van der Horst et al. 2009) , the observations at 820 MHz 
sampled more than 97% of the luminosity distribution of 
the population of known ordinary pulsars with rotational 
period longer than 100 ms, as derived from the ATNF 
pulsar cataloguer^ 



5. PLATEAU DE BURE MM OBSERVATIONS 

SGR 0418+5729 was observed with the Plateau de 
Bure Interferometer (PdBI) at 1.8 mm (166.50 GHz) in 
the D configuration between June and July 2011 (June 
27, and July 09, 10, 15, and 16). This configuration pro- 
vides baselines between 22.1 and 95.6 m. The phase cen- 
ter of the observations was 04:18:33.867, +57:32:22.910. 
The dominant track was July 15 (8h-track and excel- 
lent weather conditions). The system temperatures were 
typically in the 150 to 200 K range, and the averaged 
atmospheric precipitable water vapor was 2 mm. The 
gain calibration was performed observing the quasars 
B0552+398 and J0512+294. After calibration, the phase 
rms was 20-60°. The bandpass calibrator used was 
B0851+202. The adopted flux density for the flux cal- 
ibrator 3C273 was 16.57 Jy. Calibration and imaging 
were performed using the standard procedures in the 
CLIC and MAPPING packages of the GILDA^f] soft- 
ware. The resulting final map, obtained combining all 
the data, yields a synthesized beam size of 3'.'95 x 3'/16 
with a position angle of PA = 0.0°. The rms noise 
achieved using the full 3.6 GHz provided by the WideX 
correlator is 60 /iJy beam -1 . The primary beam of the 
PdBI at 166.50 GHz is 30 / /3. 

20 http:/ /www. gb. nrao.edu/gbtprops/man/GBTpg.pdf 

21 http:/ /www. atnf.csiro.au/research/pulsar/psrcat/ 

22 GILDAS data reduction p ackage is available at 
|http: / /www.iram. fr/IRAMFR/GILDAS 
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Fig. 7. — Left panel: Pulse profile of the first (blue) and last (black) XMM-Newton observations in the 0.5-10 keV energy band. Right 
panel: phase-resolved unfolded spectra for the first XMM-Newton observation. The spectra are relative to phases: 0.0-0.4 (black), 0.4—0.6 
(red), and 0.6-1 (green). The phase ranges are relative to the blue pulse profile in the left panel. 



We did not detect continuum emission within the 
PdBI primary beam towards SGR 0418+5729 , and ob- 
tained an upper limit of 0.24 mJy beam -1 at a 4<t 
level (see Figure[8l. The only detected source is at 
RA=04:18:30.0777"T>ec=57:32:52.00, which corresponds 
to an offset of (30.5", 29.1") with respect to the phase 
center, or a total offset of 42.1". The flux density of 
this millimeter source is 0.34±0.06 mJy (from a Gaus- 
sian fit in the uv plane and without correcting for the 
primary beam response). In addition, we looked for pos- 
sible "pulses" of emission at 1.8 mm. In order to do 
that, we checked the calibrated amplitude vs time for 
the longest track (July 15th). By averaging the visibili- 
ties in intervals of 1 minute, we found no hints of variable 
emission at an upper limit of roughly ~ 10 mJy. 

We also searched the NRAO VLA Sky Survey at 21 
cm and found no source within 15' of the millimeter 
source (Condon, et al. 1998; limiting brightness: 2.0 
mJy beam -1 ). 

6. WILLIAM HERSCHEL TELESCOPE OPTICAL 
OBSERVATIONS 

We acquired four 300s r-band images of the field con- 
taining SGR 0418+5729 on 2009 August 16, using the 
ACAM imager mounted at the 4.2 m William Herschel 
Telescope on La Palma. The average seeing was 1" and 
airmass 1.26 . Observations of a nearby field contain- 
ing SDSS calibrated stars were obtained for the absolute 
photometric calibration, while astrometry was performed 
against 2MASS sources, resulting in an accuracy of —0.1" 
on both RA and Dec. 

No source was detected within the 95% confidence 
down to a limit of r = 24. PSF photometry re- 
veals that the nearest object is detected at a magnitude 
r = 22.7 + 0.1 and center coordinates RA=04:18:34.0, 
Dec=57:32:23.5. This source is consistent with the near- 
infrared source reported by Wachter et al. (2009), but 
its distance from the SGR0418+5729 position (-1.4"), 
makes the association with the magnetar rather unlikely. 

7. FERMI-LAT GAMMA-RAY OBSERVATIONS 



We used data from the Large Area Telescope (LAT) 
onboard Fermi (Atwood et al. 2009) from 2008 Au- 
gust 4 until 2012 October 24. The Fermi Science 
Tools SC09-28-00 package is used to analyze the data. 
We selected events from the "Source" class of the 
"P7 . 6_P130_BASE" data version within a circular region 
of interest (ROI) of 10° radius centered on the position 
of SGR 0418+5729,, and in the energy range 100 MeV- 
100 GeV. The good time intervals are defined so that 
the ROI does not fall below the gamma-ray-bright Earth 
limb (defined at 100° from the Zenith angle), and the 
source is always inside the LAT field of view, namely 
in a cone angle of 66°. The "P7S0URCE_V6" instrument 
response functions (IRFs) are applied in the analysis. 

The likelihood analysis of SGR 0418+5729 was per- 
formed by means of the binned maximum-likelihood 
method (Mattox et al. 1996), using the official tool 
gtlike released by the Fermi-LAT collaboration. The 
spectral-spatial model created for the likelihood analy- 
sis includes the Galactic, and the Isotropic diffuse emis- 
sion models, as well as all the 2FGL sources within a 
radius of 15° from SGR 0418+5729. Since there is no 
2FGL source that is positionally associated to the mag- 
netar, we added in the spectra-spatial model a point-like 
source modeled with a simple power-law with the coor- 
dinates of SGR 0418+5729. The 2FGL sources within 3° 
of SGR 0418+5729 (3 sources) are modeled with the flux 
parameter allowed to vary, while the others 34 sources 
had all their parameters fixed to the value from the 2FGL 
catalog. Figure|8] shows the diffuse subtracted TS map 
of the 7° x 7° region centered on SGR 0418+5729. It was 
obtained associating to each pixel (of size 0.1° x 0.1°) the 
TS value calculated assuming a point-like testing source 
in its center. Diffuse subtracted TS map means that the 
spectral-spatial model for the null hypothesis includes 
only the Galactic and Isotropic emission models, so that 
the point-like sources should be visible in the map. 

As no significant gamma-ray counterpart to 
SGR 0418+5729 is identified, 95% flux upper limit 
is derived using the Bayesian method developed by 
Abdo et al. (2010). The 95% flux upper limit for 
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Fig. 8. — Left panel: Contours of the 1.8 mm Plateau de Bure emission of the field of SGR 0418+5729. Contours are -4, -2 (dashed grey), 
2, and 4 (white) times the rms noise of the map, 60 /^Jybeam" 1 , and they are over-plotted on the Wide-field Infrared Survey Explorer 
image at 4.6 fim. The star symbol indicates the position of SGR 0418+5729, and its size corresponds to its positional uncertainty (~ 1.2" 
in diameter). The synthesized beam, of 3.95" X 3.16", at P. A. = deg, is shown in the bottom right corner. The square indicates the field 
of view of the r-map acquired with the william Herschel Telescope (central panel). Central panel: William Herschel Telescope r-band field 
of SGR 0418+5729. Right panel: Fermi-LAT (0.1-100 GeV). Diffuse subtracted TS map of the 7° X 7° sky region centered on the magnetar 
position. The map is calculated for E > 300 MeV. The 2FGL sources are labeled in green, while SGR 0418+5729 in magenta. 



E > 100 MeV is F < 1.3 x 10~ 8 photons cm~ 2 s" 1 , 
including systematics. 

The non detection of SGR 0418+5729 at energies > 
100 MeV is not surprising, given a similar non detection 
of all other known magnetars (Abdo et al. 2010). 

8. DISCUSSION 

We have presented here a detailed X-ray study of the 
outburst of the low magnetic field soft gamma repeater 
SGR 0418+5729. The long term monitoring we per- 
formed over 1200 days allowed us to measure the period 
derivative of this pulsar (P = 4(1) x 10 _15 ss _1 ) with a 
3.5cr significance (Figure[T]). This yields an estimate of 
its dipolar magnetic field of Ba ip ~6x 10 12 G , and con- 
firms this object as the magnetar with the lowest dipolar 
magnetic field ever discovered. 

Assuming that SGR 0418+5729 attained its quiescent 
state in the last few observations, the X-ray quiescent 
emission appears dominated by a very small spot at 
kT r*> 0.3 keV and of radius ~0.16 km (assuming a 
2 kpc distance) , which corresponds to a cap of semi- 
aperture ~ l°-2°, similar to what observed in old ra- 
dio pulsars. However, in the present case the rota- 
tion power (E sa 10 29 ergs _1 ) is about two orders of 
magnitude smaller than the observed X-ray luminosity 
(« 10 31 ergs _1 ), thus indicating a different origin for the 
emission, most likely magnetic. Actually, it is quite likely 
that most of the surface is at a much lower temperature 
and is therefore invisible at energies between 0.5-10 keV. 
This implies that the quiescent bolomctric flux may be 
severely underestimated (see also below). 

The study of the spectral evolution during the outburst 
shows the presence of a non-thermal component (proba- 
bly magnetospherical) at the beginning of the outburst, 
which fades away after a few hundreds days. On the 
other hand, the temperature of the small region respon- 
sible for the surface anisotropy fades from 0.9 to 0.3 keV 



in a few year timescale (see Table 1 and Figure[3]) . 

The pulse profile evolution during the outburst decay 
shows some interesting features. As Figure[5]shows, there 
is an overall trend towards a simplification of the pulse, 
which starts with a complex, three-peaked shape and 
ends with a fairly sinusoidal pattern. Furthermore, the 
study of the pulse profiles as a function of energy (Fig- 
ure]^) in the first outburst stages, shows a great variabil- 
ity too. A large dip in an otherwise rather sinusoidal 
profile is observed at energies between 1-4 keV. 

The large pulsed fraction of 40-70% (Figure II) , the ev- 
idence of nearly phase aligned spots responsible for the 
0.9 keV thermal emission in the early outburst phases, 
as well as the 0.3 keV emission at late times, disfavor the 
presence of two spots at different temperatures, while fa- 
voring the presence of a single tiny spot cooling down 
(from 0.9 to 0.3 keV) and reducing its size (from 0.21 to 
0.16 km) during the return to quiescence. This means 
that the multi-peaked pulse profile is probably due to 
anisotropics in the magnetospheric electrons distribution 
(on top of a non- isotropic surface thermal emission). 

8.1. SGR 0418+5729 as an evolved magnetar 



In Turolla et al. (2011 1 it was shown that the rotational 
properties of SGR0418+5729 can be reproduced if the 
source is an aged magnetar, which experienced substan- 
tial field decay but still retains a strong enough internal 
toroidal field. The most updated magneto-thermal evo- 
lutionary models discussed in (Vigano et al. 2013; but 
see also Pons et al. 2009, and Aguilera et al. 2008), 
confirm this scenario. The evolution of an initial dipolar 
magnetic field of B® ip ~ 1.5 x 10 14 G (surface value at 

the pole) correctly provides the observed P and P at an 
age of ~ 550 kyr, which is probably the real age of this 
source. 

Although different combinations of the initial compo- 
nents of the magnetic field are possible, in all the mod- 
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Fig. 9. — Magneto-thermal evolution of a neutron star with an initial poloidal field of B^i p = 1.5 X 10 14 G: period (left-top), period 
derivative (right-top), evolution in the P-P diagram (left-bottom), and bolometric thermal luminosity (right-bottom). The gray band 
corresponds to the uncertainty of the angle-dependent spin-down formula. 



els the magnetic field must have been large in the past 
(> io 14 G) to explain at the same time the long spin 
period, the bright X-ray emission at this old age, and 
the flaring activity of the source. The characteristic age 
overestimates the real one by almost two orders of magni- 
tude. In Figure[9j we show the evolution of period, period 
derivative, the source track in the P-P diagram, and the 
bolometric thermal luminosity. In this scenario we esti- 
mate that SGR0418+5729's mean surface temperature 
should be now of ^0.05 keV, unfortunately undetectable 
by current X-ray observations (which are observing only 
a hot tiny region on the star surface). 

For the evolution of the timing properties, we assume 
the magneto-dipole braking formula given by Spitkovsky 

(2006): IQQ « B ^" 4 (1 + sin* 2 ), where R is the NS 
radius, x is the angle between the rotational and the 
magnetic axis, c is the speed of light, = 2ir/P is the 
angular velocity, and I is the moment of inertia of the 
star. 

An alternative possibility is that the neutron star was 
born with an external magnetic field close to the present 
one, but its large core poloidal field slowly diffuses out 



i.e. by ambipolar diffusion (Soni 2012). Although from 
the timing properties alone it is hard to discriminate be- 
tween a hidden strong crustal magnetic field, a hidden 
strong core field or an intrinsically low-B neutron star, 
the magnitude of the X-ray luminosity, and the spectral 
properties and light curves may be used to distinguish 
the different scenarios. 

In particular, the low magnetic field scenario cannot 
explain the high luminosity, large pulsed fraction, and 
the flaring activity of the source. As a matter of fact, 
if there is little field decay and the real age corresponds 
to the characteristic age (which in this scenario would 
be needed to reach the present period of 9 s), no existing 
non-magnetic cooling model can account for an X-ray 
luminosity of w 10 31 erg s _1 at a characteristic age of 
« 35 Myrs. 

The scenario in which a large core field diffuses out 
has the same problem: if the real age of the star is simi- 
lar to its old characteristic age, no cooling model in the 
literature predicts such high quiescent luminosity. In 
addition, while the timescales used in Soni (2012) are 
correct for normal, non-superfluid nuclear matter, recent 
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work (Glampedakis, Jones & Samuelsson 2011, MNRAS) 
shows that in the presence of superfluidity in the neutron 
star core, the timescales for ambipolar diffusion are many 
orders of magnitude longer, and therefore ambipolar dif- 
fusion does not play any role during the active age of the 
star. 

8.2. SGR 0418+ 5729 outburst rate 

An important question is whether a relatively low dipo- 
lar field is consistent with the star-quake model in which 
the primary cause of the outburst is an internal depo- 
sition of energy following a crust fracture. It is often 
overlooked that the magnetic stress needed to break the 
crust is strongly dependent on density (it is much easier 
to break the outer crust than the inner crust) and that 
the crust thick ness grows as the temperature drops with 
age. In Figure [T0| we show an estimate of the minimum 
magnetic field variation required t o induce a fracture. 
As assumed in Perna & Pons (2011 ), this estimate is ob- 
tained assuming that the crust moves through a series of 
equilibrium states in which its elastic stress balances the 
(time-dependent) magnetic stress. The deviation of the 
magnetic field with respect to the last unstressed con- 
figuration (SB C ) may be large enough to break the crust 
when 

sb c « (K") 1/2 (i) 

where cr™ ax is the maximum stress that a neutron star 
crust can sustain (Chu et al. 2010). For young, rela- 
tively hot magnetars (crustal temperatures of 5 x 10 8 K), 
only the inner crust is solid, and strong field variations 
SB C > 10 14 G are required to fracture the crust. How- 
ever, for old, cold neutron stars (crustal temperatures of 
10 8 K), the solid crust extends down to 10 8 g/cm 3 , and 
is much easier to break, even with variations of the mag- 
netic field of the order of SB C > 10 12 G. Note also that 
fractures close to the inner crust are much more energetic 
(because of both the higher available elastic energy and 
larger volume involved) than fractures in the low density 
region. In the first case, one can reach up to 10 44 erg, 
while in the second case, events of ~ 10 41 erg are ex- 
pected. A rough pred i ction of the expected outburst rate 
( |Perna & Pons||2011| |Pons fc Perna||2011[ ) for the solu- 
tion model mentioned above gives < 10 -3 star-quakes/yr 
for an object as SGR 0418+5729. Assuming that there 
are about 10 4 neutron stars in the Galaxy with similar 
age, and that a (very approximatively) 10% of them are 
born as magnetars, a naive extrapolation of this event 
rate to the whole neutron star population leads to an 
expected low-B magnetar outburst rate of < 1 per year. 
Therefore, we expect that more and more objects of this 
class will be discovered in the upcoming years (as e.g. 
Swift 1822.3-1606; Rea et al. 2012; Scholz et al. 2012). 

8.3. SGR0418+ 5729 outburst decay 

By modeling the flux evolution in time we teste d if the 
crustal cooling model presented in |Pons &: Rea ( 2012[ ) 
can fit the flux decay of SGR 0418+5729, on the wave 
of what done for Swift J1822. 3-1606 (Rea et al. 2012). 
We assume a dipolar field of 6 x 10 12 G (equatorial), an 
int erna l toroidal field of 10 14 G (at maximum) as inferred 
in §8.1 and an average surface temperature of 0.05 keV, 
which is the temperature we expect for the surface of 



such an old magnetar (note that in the 0.5-10 keV band 
we are only seeing a tiny hot spot). The best modeling 



was found by injecting 2.5 x 10 26 erg cm" 3 in a thin layer 
in the outer crust between 4.5 x 10 9 and 10 10 g cm -3 , and 
in the region contained within a cone with axis in the di- 
rection of the magnetic pole, and aperture a« 0.4 rad, 
for a total energy deposition of 2.5 x 10 41 erg (compatible 
with typical magnetar outbursts; Pons & Rea 2012). The 
evolution of the bolometric, and of the 0.5-10 keV flux is 



shown in Figure 10 (solid and dashed lines, respectively). 
We have shown with red squares how the observed flux 
decay would appear when adding the contribution of a 
blackbody component at 0.05 keV, mimicking the entire 
neutron star surface. It is clear that crustal cooling can 
easily explain the decay only if this further component is 
taken into account. In particular, the solid line is fitting 
the red points because the entire neutron star surface is 
taken into account, while it is not in the observed black 
data, which in fact cannot be reproduced by the dashed 
line. This is indicative of the difficulty of comparing the- 
oretical cooling curves with data obtained in a certain 
energy band. We also note that no theoretical model 
predicts a surface temperature as high as 0.3 keV on a 
timescale of years, unless a continuous energy release is 
assumed (e.g. by long-lived internal currents). We fi- 
nally mention that all the previous considerations are 
based on the (implicit) assumption that the blackbody 
temperature is a measure of the physical temperature of 
the emitting region. If this turned out not to be the case, 
e.g. because the spectrum is thermal but not Planck- 
ian so that a color correction is required, the physical 
surface temperature may be smaller than the measured 
blackbody temperature. 

An alternative model to the crustal cooling scenario 
consists in the presence of currents flowing into the mag- 
netosphere through a gradually shrinking magnetic bun- 
dle heating the neutron star surface from the top. In 
particular, the deepest available XMM-Newton observa- 
tion of SGR 0418+5729, performed two months after the 
outburst onset, reveled that the 0.5-10 keV spectrum of 
SGR 0418+5729 is best reproduced by a blackbody com- 
ponent plus an additional non-thermal component, or by 
a resonant cyclotron scattering model. This suggests the 
presence of twisted magnetic field lines, at least in the 
first outburst stages. Furthermore, the limited spatial 
extent of the heated region (< 1km) is suggestive of a 
scenario in which the twist is confined within a small part 
of the ma gnetosphere, a thin current-carrying bundle, or 
j-bundle ( |Beloborodov 20091. As the j-bundle untwists 
during the outburst decay, the spectrum becomes more 
and more blackbody-like, as indeed observed. 

Resonant cyclotron scattering from a thin, decaying j- 
bundle appears also capable of explaining (qualitatively) 
the spectrum and its evolution during the first outburst 
stages, but whether it can explain the double-peaked 
pulse profiles is unclear. However, this scenario has 
some further difficulties: 1) the total luminosity pro- 
duced by currents in the bundle is, for such a low-i? 
and a small thermal spot, well below the one observed at 
early times, at least if th e spot is at the polar cap (se e 
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Beloborodov 2009 



2011) 



and also Turolla et al 
2) the timescale lor the twist decay is much shorter 
(< 1 yr) than what implied by the long outburst of 
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Fig. 10. — Left panel: Minimum variation of the magnetic field required to break the crust by magnetic stresses as a function of density. 
The vertical dashed lines delimit the transition from solid to liquid for three different temperatures. Right panel: Outburst modeling 
following Pons & Rea (2012). Black data are the 0.5-10 keV unabsorbed flux, while red squares are the bolometric unabsorbcd flux with 
the addition of the flux of a thermal component at kT=0.05keV from the entire neutron star. Solid and dashed lines refer to the outburst 
model for the bolometric and 0.5-10 keV thermal flux, repectively (see text for details) 

SGR 0418+5729, and 3) an approximate relation between 
the emitting area and the luminosity exists (A ~ L 2 ; Be- 
loborodov 2009) if most of the luminosity is produced 
by current dissipation, but SGR 0418+5729 data show a 
somewhat natter dependence when the first stages of the f 
outburst are fitted (see Figure 11). | 
In summary, the crustal cooling model, when including J 
also a possible hidden contribution from the entire neu- t 
tron star surface, appears favourable in explaining the « 
outburst decay of SGR 0418+5729. However, it is likely 
that a combination of crustal cooling and magnetospheric 
untwisting bundle can be operating at the first stages of 
the outburst. 




8.4. Constraints on the presence of a fossil disk 
surrounding SGR 0418+5729 

A fallback disk around SGR 0418+5729 was suggested 
by Alpar et al (2011) as a way to aid the spin-down of the 
pulsar and explain the 9 s periodicity of t his s ource. As 
an alternative, our results show (see Figure 8.1 ) that both 
the thermal and the timing properties of this source can 
be reproduced for a pulsar age of ~ 550 kyr by properly 
accounting for magnetic field evolution and dissipation, 
which also imply that the neutron star was born with a 
much hig her dipolar field than the one measured today 
(see 5 8.1). Hence, in principle, the timing properties of 
this source would not necessarily require an additional 
spin down torque by a disk. However, given the sug- 
gestion that fallback disks around isolated neutron stars 
might be common (Michel 1988; Chevalier 1989; Lin et 
al. 1991), it is worthwhile to use the current multi-band 
upper limits to set constraints on the presence of a fall- 
back disk around SGR 0418+5729. 

Since the pulsar is currently spinning down, any 
disk-magnetosphere interaction must probably occur in 
the propeller regime, with the pulsar transferring an- 
gular momentum to the disk. For this condition to 
be satisfied, the inner boundary of the disk, located 
at about the magnetospheric radius R m = 2.5 x 
10 8 [M/(10 16 gs- 1 )]- 2 / 7 (M NS /M )- 1 / 7 [ J B/(lO 12 G)] 4 / 7 , 
must be equal to or larger than the corotation radius 
Rco = (GMns) 1 / 3 ^ 2 / 3 . The strongest constraint on the 




Fig. 11. — X-ray luminosity evolution as a function of the black- 
body emitting area and temperature (see also Table 1; we assume 
a 2kpc distance). 

disk emission is obtained when the inner radius of the 
disk obtains its minimum value, i.e. R ln — R m = R co . 
For the outer radius, we assume R out — 10 14 cm. We 
found that larger values do not result in appreciably 
larger emission at the frequency of interest, and hence 
this value allows to set the tightest constraint on the 
disk emission. 

With the inner and outer disk radii fixed as discussed 
above, the emission spectra from the fallback disk is com- 
puted using the model of Perna et al. (2000). The disk is 
assumed to be optically thick and geometrically thin, and 
the anisotropy in the X-ray luminosity from the source 
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(which irradiates the disk) is neglected, since it is found 
to be of second order (Perna & Hernquist 2000). The 
disk is assumed to be still "active", i.e. viscously accret- 
ing (see Menou et al. 2001). The disk emission is the 
result of both viscous dissipation and re-radiation of the 
pulsar X-ray luminosity. In order for the magnctosphcric 
radius not to exceed the corotation radius, the accretion 
rate must be limited to M ^ 10 15 g s _1 . With this value, 
the disk luminosity in the mm band is dominated by re- 
processing of the pulsar X-ray luminosity. At 166 GHz, 
the predicted flux is about 0.01 mJy for a face-on disk, 
below the measured limit of 0.24 mJy. Hence the pres- 
ence of a fossil disk cannot be ruled out by the current 
mm measurements. Even adding the contribution from 
the whole surface of the star by a putative thermal com- 
ponent at 0.05 keV, would bring the predicted mm flux 
just around the measured flux limit for a face-on disk. 

The field around SGR 0418+5729 was also observed 
with the Grantecan and Hubble telescopes (Esposito et 
al. 2010; Durant et al. 2011). In particular, the 
latter observations were performed in two wide filters, 
the optical, with a pivot wavelength of 5921 A, and 
in the NIR, with pivot wavelength of 11534 A. The 
source was not detected down to the flux of fo < 
2.3 x 10~ 31 erg s" 1 cm~ 2 Hz" 1 and f NIR < 4.4 x 
10~ 31 erg s -1 cm -2 Hz -1 , respectively. We found that 
this optical limit (nor the Grantecan or WHT limit) is 
not sufficiently constraining for a disk with the proper- 
ties described above (the predicted emission for a face-on 
disk is about a factor of four below the limit). On the 
other hand, in the NIR, the observational limit is al- 
ready able to rule out a face-on disk, which would yield 
an emission about twice larger than the measured flux 
limit. However, a disk inclined with respect to the ob- 
server by an angle cos 9 £ 0.5 would still be allowed by 
the observations (although falling short in explaining the 
X-ray bursts of this object). 

8.5. Conclusions 

At the time of writing, in the ATNF pulsar cata- 
logue (Manchester et al. 2005) 138 isolated radio pul- 
sars have a dipolar magnetic field larger than that in- 
ferred for SGR 0418+5729. Our results imply that some 
of these objects might hide a strong toroidal component 



of the internal field, not measurable via the pulsar tim- 
ing properties. A hint for such strong fields might be 
a high surface temperature, hotter than what would be 
predicted by standard cooling models at the pulsar age. 
However, only a few of those pulsars have had dedicated 
X-ray observations, and the shallow surveys do not suf- 
fice to detect such emission (expected to be as luminous 
as L x ~ 10 31 crgs _1 ). Furthermore, our calculation of 
the outburst rate of a low magnetic field magnetar also 
suggests that roughly once a year a quiet neutron star 
might turn on with magnetar-like activity. 

On the other hand, if indeed a large number of neutron 
stars is hiding a strong magnetic field component, there 
would be important consequences also for other branches 
of astrophysics. In particular, it would imply that super- 
nova explosions should generally produce strong mag- 
netic fields, and that most massive stars are either pro- 
ducing fast rotating cores during the explosion to acti- 
vate a dynamo, or are strongly magnetized themselves. 
Furthermore, in this scenario a non-negligible fraction of 
gamma-ray bursts might be due to the formation of mag- 
netars, and the gravitational wave background produced 
by magnetar births should then be larger than predicted 
so far (important for future instruments as Advanced- 
LIGO). 
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